Little-Parks oscillations near a persistent current loop 
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We investigate the Little-Parks oscillations caused by a persistent current loop set on the top 
edge of a mesoscopic superconducting thin-walled cylinder with a finite height. For a short cylinder 
the Little-Parks oscillations are approximately the same ones of the standard effect, as there is 
only one magnetic fiux piercing the cylinder. For a tall cylinder the inhomogeneity of the magnetic 
field makes different magnetic fluxes pierce the cylinder at distinct heights and we show here that 
this produces two distinct Little- Parks oscillatory regimes according to the persistent current loop. 
We show that these two regimes, and also the transition between them, are observable in current 
measurements done in the superconducting cylinder. The two regimes stem from different behavior 
along the height, as seen in the order parameter, numerically obtained from the Ginzburg-Landau 
theory through the finite element method. 
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I. INTRODUCTION 

In the last decades we have witnessed a huge progress 
in the microfabrication of systems and measurement 
techniques that have allowed the study of mesoscopic su- 
perconducting structures^. The critical temperature of a 
superconducting ring oscillates according to the applied 
external magnetic field, as shown by Little and Parks 
(LP) in 1962^. The LP effect can be regarded as a fore- 
runner of recent developments^, because of its eminently 
mesoscopic nature. It predicts that for a thin-walled 
cylinder of radius R, the critical temperature varies as 
ATc = {h^/8mR'^){^/^o - ?^)^ where m is the Cop- 
per pair mass, $o is the fundamental flux, $ = HttR^ 
is the magnetic flux trapped inside the ring, and n de- 
scribes a quantum number. A temperature variation de- 
tectable within experimental range, say ATc ^ 10^^ K, 
means that the radius must be in the mesoscopic do- 
main, R ~ 1.0/^m. The LP oscillations are detected 
by measuring two consecutive temperature maxima, each 
occurring for ATc = 0, when the total magnetic flux $, 
which is the sum of the external magnetic flux piercing 
a given surface and the magnetic flux produced by the 
circulation of the screening supercurrents along the curve 
bounding this surface, add up to n$o- Since its initial 
proposal the LP effect has been measured by several dif- 
ferent techniques and several distinct systems, such as a 
perforated disk with varying hole size*, a patterned mi- 
crostructure of an oxide superconductor—, a single meso- 
scopic Al ring (SQUID measurement of the susceptibil- 
ity )S, an array of Al loops (specific heat measurement )i, 
and a single YBCO sub-micron ringS^. The LP effect can 
be used to demonstrate many features of the supercon- 
ducting state, such as the interaction among two super- 
conducting order parameter o^do and the study of super- 
conductor/ferromagnet hybrida^i. The LP effect is well 
described by the Ginzburg-Landau theory in most of the 
situations^"—. The LP can also exist through fluctua- 



tions in the normal phase for very small rings with effec- 
tive radius R/^ < 0.6 under a constant magnetic ficld^^. 

In this paper we investigate the LP effect in a thin- 
walled cylinder with a flnite height produced by the in- 
homogeneous magnetic held of a persistent current loop, 
put on its top. This external magnetic held is that of a 
magnetic moment pointing along the thin-walled cylin- 
der major axis. In the standard LP effect the height of 
the thin-walled cylinder plays no role because the ap- 
plied external field is constant and oriented along the 
major axis. However in the presence situation the height 
plays a major role^^. We report here new features due to 
the inhomogeneity distribution of the supercurrent in the 
cylinder, which in its top edge feels a magnetic flux more 
intense than in the bottom. Our theoretical framework 
is developed in a current bias regime, which means that 
the persistent current circulating in the loop is varied and 
the response of the thin-walled cylinder observed. 

The persistent current loop of a (non-superconducting) 
metallic mesoscopic ring has been observed since decades 
ago^^ and in principle can be used for an experimental 
realization of the present proposal. The persistent cur- 
rent loop created by a superconducting mesoscopic ring 
is another way to realize the present system. In this case 
the system has two superconductors, for instance, Nb 
and Al, the former used for the persistent current loop 
and the latter for the the thin- walled cylinder. Since 
the Nb critical temperature (Tc = 9.25 K) is 7.8 times 
larger than that of Al (T^ = 1.18 K)i^, the LP oscilla- 
tions observed in the thin- walled cylinder are near to the 
Al transition to the normal state, which happens to fall 
in the very low temperature domain of the Nb current 
loop. Therefore the Nb ring generates a stable steady 
magnetic whose value is swept continuously from zero 
to its maximum value by varying the persistent current 
to the maximum Nb critical current. Let a and R be 
the radii of the external current loop and of the thin- 
walled cylinder, respectively, as shown in Fig. [TJ We 
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are mostly interested here in the effects brought by the 
inhomogeneity of the field, and for this reason will con- 
centrate in the case R > a. The magnetic moment of 
the persistent current loop is given hy ^ — Iioop t^o? but, 
notice that the thin-walled cylinder also defines a mag- 
netic moment scale, given by /iq — $0'?o/27r, where 
is the zero temperature coherence length of the super- 
conductor that makes the Al thin- walled cylinder. We 
study the LP osciUations by sweeping the ratio /x//xo to 
its maximum value ^c/mo- Thus it is of fundamental im- 
portance to determine the maximum attained magnetic 
moment in the current loop, Hc = Ic tto^, Ic being the Nb 
critical current. According to Sebastien Michotte et alJ^ 
for T= 1.5 K the Nb critical current is Ic=2000 ^A for a 
wire of width 75 nm. Therefore a ring with radius a=0.5 
/im achieves the maximum magnetic moment /ic=1.5 x 
10^^^ A.m^. To determine the magnetic moment scale 
/io of the thin-walled cylinder it suffices to know the zero 
temperature coherence length of Al. Nanoscopic Al rings 
have been extensively studied in the literature, such as 
by Florian R. Ong et alii and Hendrik Bluhm et al.^. 
Data fit in both cases leads to the values of equal 
to 150 and 70 nm, respectively. Here we shall consider 
^0=100 nm which gives ^o=3.2 x 10~^^. In summary we 
find that maximum magnetic moment ratio is /ic/Aio=3 
X 10^. Therefore the current circulating in the external 
loop can be swept up to six orders of magnitude in our 
suggested system. Notice that our theoretical calcula- 
tions do not rely in this specific realization which is only 
a suggested experimental set up for it. 

To describe the present system two other lengths must 
be considered, namely, the distance from the current loop 
to the center of the thin- walled cylinder, zi and the height 




FIG. 1. (Color online) This figure depicts the persistent cur- 
rent loop (radius a = S.O^o) on top of a thin-walled super- 
conducting rod (radius R — lO.O^o) put 0.5^0 away from its 
edge. Two heights are treated here, corresponding to a short 
{Zo = 2.05o) and a tall (10.05o)cylinder. 



of the thin- walled cylinder, = 2zo, zq being the dis- 
tance from the edge to its center. Throughout this paper 
the distance between the current loop and the top edge 
of the thin- walled cylinder is kept fixed and equal to half 
coherence length, thus zi — 0.5^o + zq- Our analysis is 
restricted to the thin-walled cylinder, which means that 
its width is smaller than the coherence length, w < 
and therefore only two coordinates are needed to deter- 
mine a point in the cylinder, namely, {p — R, Lp, \z\ < zq) 
in cylindrical coordinates with the origin at the center 
of the thin- walled cylinder. The order parameter inside 
the thin-walled cylinder is expressed by its phase and 
amplitude /, such that 5' — f{z) exp {i(j){ip,z)). The well 
known condition of a single value order parameter implies 
that '^{ip -f 27r) — ^'((ys) and this means that the order 
parameter must be the same along the z-axis: (j) = Lip, 
where L is an integer to be associated with the angular 
momentum of the vortex state trapped in the thin- walled 
cylinder. Therefore the order parameter is assumed to 
have the general form 'i/{z,(p) — f (z) exp (iLip). The 
dipolar magnetic field falls with the inverse of the cube 
of the distance to the current loop, thus being strong in 
the top and weak at the bottom and this will make the 
amplitude strongly z dependent. Although the current in 
the external loop always makes the thin-walled cylinder 
feel a dipolar magnetic field, its inhomogeneity is only 
noticeable in case the thin-walled cylinder is sufficiently 
tall. For this reason we numerically study here two ex- 
amples of thin- walled cylinders here, a short (Zq = 2.0^o) 
and a tall {Zq — lO.O^o) one. The former basically yields 
standard LP oscillations and the latter comprises novel 
phenomena because the top and the bottom edges of the 
thin- walled cylinder are pierced by very distinct magnetic 
fiuxes. 

We solve the nonlinear GL equation for the thin- walled 




FIG. 2. (Color online) The LP temperature Td/J.) versus 
the magnetic moment fi is shown here for the short cylin- 
der (height Zo = 2.0^0 )• The L lines, ranging from to 20, 
are obtained from the linear theory fEa. (|13|) ) and describe 
angular momenta trapped in the cylinder. 
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cylinder through the finite element method (COMSOL 
and MATLAB softwares) and check that our numerical 
procedure is correct in the short cylinder limit, since in 
this case, we also solve the theory analytically. This gives 
confidence that our numerical results for the thin- walled 
cylinder are reliable. 

The paper is organized as follows. In section we 
review the basic Ginzburh-Landau formalism applied to 
the thin- walled cylinder in present of a persistent current 
loop. In section Hill we explore the limit that the ampli- 
tude / can be considered as constatn along the height to 
obtain important theorems relating the zero current to 
the minimum of the free energy. In section ITVl we numer- 
ically study the short and tall cylinders and find that the 
latter displays two distinct regimes according to the the 
value of the persistent current. Finally in section fVl we 
reach some conclusion about the results obtained here. 



II. THEORETICAL BACKGROUND 



where a{T) = ao ^y- — 1 j , therefore a{T) < to sustain 
the superconducting state, and /3 is a constant. The vol- 
ume of the thin-walled cylinder is 17 = 2'kRwZo and the 
dimensionless volume element is dv = (Px/il. From the 
variational principle one obtains the Ginzburg-Landau 
equation that determines the order parameter which min- 
imizes this free energy: 




= -a(T)^' -t- /Jl^'p^' (2) 



The persistent current loop has its center located at co- 
ordinates (0, 0, zi), its magnetic moment is /i = IioopT^a^, 
hoop being the circulating current. The vector poten- 
tial describing its field in space is given by A{p, z) = 



The Ginzburg-Landau theory describes the supercon- 
ductor near to its critical temperature Tc- Because the 
London penetration length is taken to be much larger 
than the hollow cylinder height, the shielding of the mag- 
netic field is safely ignored in our calculations. Hence the 
free energy density is simply that of the sum of terms in- 
volving the order parameter: 
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FIG. 3. (Color online) The LP temperature Tc{p) versus the 
magnetic moment p is shown here for the tall cylinder (height 
Zo ~ lO.O^o). The L lines, ranging from to 41, are obtained 
from the linear theory (Eq.([T3])) and describe angular mo- 
menta states trapped in the cylinder. A few selected curves 
are shown as thick lines and their properties discussed in the 
text. 



FIG. 4. (Color online) The free energy of the short cylin- 
der (height Zo = 2.0^o) is shown here for two temperatures 
T = O.STc and T = 0.9Tc. Free energy L lines ranging from 
to 20 are shown versus the magnetic moment p of the per- 
sistent current loop. A dashed vertical line indicates that the 
position of the minima is approximately the same for the two 
temperatures. 
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The thin-walled condition restricts the radial coordinate 
to p = R, and so, it suffices to use the expression, 



A^{R, z) — —r]{z), where 



(4) 



1 



■(1 + 1.8750p2 



+4.9219^" + 14.6631p^ + 46.7386p^ + 155.4058^^° + 
+531.8949pi2 _,_ 1 8593 iQi^u ^ q qq^2 10^^ + 

+2.3757 lOY^ + 8.6335 lOV" + ■ ■ • ) (5) 



P = 



i + (i)' + (Tr^)'' 



(6) 



A term C2nP^" of this series expansion has the coeffi- 
cient C2n growing with n, as shown above, that must be 
overcome by p^" in order to achieve convergency. For 
instance, in the present case of interest, R = lO.O^o, 
a = 5.0^0, and we take the top edge of the cylinder, 
z — zi = —0.5^0 to obtain that p = 0.3992. Preci- 
sion of the order 10"'^ is only achieve near the n — 10 
term when C20 = 8.633510^ times p^° = 1.05 10"^ gives 
C2op20 = 9.0651 10-*. 

It is worthwhile to compute the magnetic flux that 
pierces a ring resulting from the intersection of the thin- 
walled cylinder with a plane z: 



^iti,z) 



j( dl- A{R, z) 



27r/i 



27rRA^{R,z) = —^v{z) 



R 



/Uo R 



(7) 
(8) 



Eq.® shows that the ratio between the magnetic mo- 
ments, n/no, together with the geometrical parameters, 
defines the ratio between the z plane magnetic flux and 
the fundamental flux. In order to gain insight into this 
series expansion, let us analyze it in the case treated here, 
namely, the ratio between the current loop and the thin- 
walled cylinder radii is equal to one half, a/R = 0.5. 
Suppose a point at the surface of the thin-walled cylin- 
der near to its top, such that (2 — zi)/R « 0. In this 

case we have that « 0.16 and 1 -I- -I- ^P^ H ~ 

1 -I- 0.30 -I- 0.12 -I- • • • . We learn that the series terms con- 
tribute in this example, although they do not add any 
signiflcant qualitative change if instead we had consid- 
ered the point like dipole limit a — 0. Thus in this limit 
all the corrections in powers of (a/i?)^ do not exist and 
the field is that of a point like magnetic dipole, fl = fj,z 
given hy A = jl X f'/\ff' where r = Rp + (z — zi)z falls in 
the surface of the thin-walled cylinder. 

We write the order parameter in terms of its amplitude 
/ and phase cj) — Lip to obtain the following expression for 
the kinetic energy in terms of the geometry, the angular 
momentum L and the vector potential A^: 



2m J n WdzJ i?2 



L - 



<P{p,z) 



'f} (9) 



Then the free energy density becomes. 



J-^G+ / dv 
Jn 



(10) 



We also write the supercurrent density, J ; 
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in terms of the amplitude and the phase de- 
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FIG. 5. (Color online) The current circulating in the short 
cylinder (height Zq = 2.0^o) is shown here for two temper- 
atures T = 0.8Tc and T — 0.9Tc. Current L lines ranging 
from to 20 are shown versus the magnetic moment of the 
persistent current loop. A dashed vertical line indicates that 
the vanishing of the current is approximately the same for the 
two temperatures. 
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FIG. 6. (Color online) The free energy of the tall cylinder 
(height Zo = lO.O^o) is shown here for the temperature T = 
O.STc. Free energy L lines ranging from to 20 are shown 
versus the magnetic moment p of the persistent current loop. 
Some L states are depicted as thick lines to exemplify the 
three existing regimes: single well (L = 3), transition (L = 7), 
and double well (L = 12). The inset helps to visualize the 
transition regime. 
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composition using the cylindrical symmetry: 



L 



(11) 



There has been many proposed methods to calculate the 
current of the LP effec t^^'^^ . In our numerical calcula- 
tions the current circulating in the thin- walled cylinder is 
simply obtained as an integration of the current density, 
J, along the cross-section of the cylinder. According to 
the thin- walled condition w ^ S^q, and so, 



r-zO 

W I dzj ■ (f 

zO 

^° dz 



Zq a/ooT/S 



L - 



(12) 



in terms of the current parameter Iq = 
(^qTiZowao) / [mf3R) . The LP oscillations occur near 
to the critical temperature, where the order parameter 
is weak, and so, besides the non-linear, the linear 
version of Eq.Q, with the cubic term neglected, can be 
used. Then there is a Schrodinger like equation for the 
amplitude f{z), once the order parameter is taken as 
'^iz,'p) = /(2)exp (iLip): 



where the potential is given by 



MT)fiz) 



(13) 



and the boundary condition at the edges of the thin- 
walled cylinder are d/(z)/dz|(2=±z(,) — 0, which is au- 
tomatically satisfied in the limit of a constant /. Let 
us consider the ill-defined limit, but still instructive, of 
/ z independent, such that d^f{z)/dz'^ « 0. Then the 
above equation defines a condition under the simple as- 
sumption that there is a non-zero /: V{z) -\- a{T) = or 

equally, T/T, = 1 - (^o/i?)' - ^(m, ^)/*o] This de- 
fines the LP temperature, which is a function of the fiux 
that pierces the thin-walled cylinder along the z plane. 
However because $(/i, z) is z dependent, the found ra- 
tio T/Tc is in contradiction with the original argument. 
Indeed a z independent amplitude / is only possible for 
a very short thin-walled cylinder, and this entitles us to 
take z = 0, which means the medium plane. We choose 
to express the above ratio T/Tc in terms of the T = 

coherence length, = ^ /2maQ. The dimensionless 
ratio $(/i, z)/$o can be written in terms of the dimen- 
sionless length ratios, R/£,o, a/R, and {z — zi)/R, and 
the dimensionless magnetic moment ratio, /i//io where ^0 
has been previously defined. This naive approach to the 
LP oscillations is improved in the next sections by treat- 
ing the z dependence of this problem. The thin-walled 
cylinder in presence of the inhomogeneous field produced 
by the persistent current loop field makes the amplitude 
dependent on its height position: f{z). In this case we 
also solve the above Schroedinger equation in order to 
determine f{z). 
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FIG. 7. (Color online) The free energy of the tall cylinder 
(height Zo = lO.O^o) is shown here for the temperature T = 
0.9Tc. Free energy L lines ranging from to 18 are shown 
versus the magnetic moment /x of the persistent current loop. 
Some L states are depicted as thick lines to exemplify the 
three existing regimes: single well (L = 3), transition (L = 5), 
and double well {L — 8). The inset helps to visualize the 
transition regime. 



FIG. 8. (Color online) The current circulating in the tall 
cylinder (height Zq = lO.O^o) is shown here for the temper- 
atures T = O.STc and T = O.QTc. Current L lines ranging 
are shown versus the magnetic moment fi of the persistent 
current loop for the L values of Figs.® and ([7]). The single 
and double well regimes reflect in current lines with different 
slopes and the two regimes are separated by a transition line, 
L = 7 for T = O.STc, and L = 5 for T = 0.9T^, respectively. 
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III. THE RING LIMIT 

A thin-wallcd cylinder such the order parameter am- 
phtude / is constant along the z direction, and therefore 
df /dz = 0, we define as a ring. The validity of the ring 
treatment restricts the height to a few coherence length 
units, which means that the ratio ^o/Co cannot be very 
large. In the ring limit all that is necessary to determine 
the order parameter and the free energy is to compute 
the integrals. 



9{^J) = 





dz 


L ~ 




J-zO 


'Z'o 




nzO 


dz 


L - 




J-zQ 




$0 -1 



and, 



(15) 
(16) 



The kinetic energy energy becomes G = 010(^)^(7/^ and 
the current / = ^^Pg'. Then, we have for the free 
energy, 

-^=K(|)'5 + «(T)]/^ + fr (17) 

whose minimization determines the order parameter and 
its minimum: 



/^ = -^K(|)^5 + «(T)] 
•^--^K(|)^5 + «(T)]^ 



(18) 



Notice that since g > and the above equation de- 
termines that ao(^)"g + a(r) < because /q > 



then there is a temperature bound given by a{T) < 

— Q!o(^) 5- The upper bond defines the LP temperature, 
which corresponds to the vanishing of the order param- 
eter, such that the superconductor reaches the normal 
state for a temperature which depends on g and the ao 
parameter, according to the above equation. Therefore 
the LP oscillations obtained from the non-linear theory 
in case of / constant are expressed in terms of the inte- 
grals g and g'. The corresponding free energy, current 
and critical temperature are defined as follows: 



Tifi, T) = - Jo 



T - r,(/i) 



-l-(|f^^(M) 



Tr 



5'(m) 



(19) 
(20) 
(21) 



The integrals g(/i) and g'(/i) can be expressed in terms 
of the following defined integrals: 



dz <i>(^, z) 

-zo Zo ^0 



and, 



zO 



dz /<i>(/^,z)\2 



-zO ^ ^0 

such that, the two integrals can be written as, 

r-zO 



(22) 
(23) 



(L) = ^h, r ^V{z), and, (24) 
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FIG. 9. (Color online) Vertical lines associated to the vanish- 
ing of the current are plotted together with the free energy 
for the short cylinder (height Zq = 2.05o) at the temperature 
T — O.STc- L states ranging from to 20 are considered here. 
The /i values equal to 9.62, 66.64, 114.04 and 190.29, labeled 
as (1), (7), (12) and (20), respectively, are depicted as red 
dashed vertical lines and are selected for further analysis in 
Fig. dm). 



FIG. 10. (Color online) Vertical lines associated to the van- 
ishing of the current are plotted together with the free energy 
for the tall cylinder (height Zo — lO.O^o) at the temperature 
T = O.STc. L states ranging from to 20 are considered 
here. The /i values equal to 12.75, 85.64, 129.02 and 264.97, 
labeled as (1), (6), (7) and (11), respectively, are depicted as 
red dashed vertical lines and are selected for further analysis 
in Fig.inj. 
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FIG. 11. (Color online) The current density J, the pair density and the potential V of the linear equation fEq. (|44|) '). 
along the cylinder height are shown here for the short cylinder (height Zo — 2.05o) at the temperature of T = O.STc. The jj, 
values, selected in Fig.®, and equal to 9.62, 66.64, 114.04 and 190.29, labeled as (1), (7), (12) and (20), are considered for this 
purpose. The persistent current loop is located at position 1.5^o. The L lines and their corresponding figures are also listed 
here. 



By taking that {^f /f ^ [t]{z) - /i] ^ > 0, we obtain 
that 

h - if > 0. (26) 

Next we show three important features of the LP oscilla- 
tions with respect to the magnetic moment fj, that can be 
expressed in terms of the above functions g{fi) and g'{fJ-), 
written in the following way: 

5(M)-i2-2L-^/i + (-^)'/2, and, (27) 

g'if,)=L-^h. (28) 
Mo 

These three three features are valid to all temperatures, 
and for this reason assume that T ^ Tc{fJ,). They are as- 
sociated to three special moments, fif,fic and fit, defined 
as such: 

1) The minimum of the free energy. According to 
Eg. dTO)) . dT/dfi = is the same condition as the 
minimum with respect to the critical temperature, 
dTc{fJ,)/dfi = 0, which boils down to dg{fi)/dfj. = 0. 



This defines the magnetic moment /i/. 

^=i^ (29) 
Mo h 

2) The vanishing of the current. According to Eq. (|2ip 
the condition /(/i) = happens for g'(/i) = 0. This 
defines the magnetic moment fic- 



Mo h 

3) The highest LP temperature. According to Eq. (j20p 
this occurs when the LP temperature becomes the 
critical one. However the condition Tc{n)/Tc = 1 
is the same as g{iJ,) — 0. This defines the magnetic 
moment fit- 
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FIG. 12. (Color online) The current density J, the pair density \'ip'^ and the potential V of the linear equation (Eq.(|44])), along 
the cylinder height are shown here for the tall cylinder (height Zo = lO.O^o) at the temperature of T = O.STc. The /i values, 
selected in Fig.((Tni), and equal to 12.75, 85.64, 129.02 and 264.97, labeled as (1), (6), (7) and (11), are considered for this 
purpose. The persistent current loop is located at position 5.5^o. The L lines and their corresponding figures are also listed 
here. 



Then we prove that the LP temperature never reaches the 
critical temperature Tc, unless for the zero height ring, 
where $0(^,2^) = *o(MjO) {h = iD- Notice that 



h 

= — > 1 
^1 



(32) 



which means that the current in the thin-walled cylinder 
always vanishes after the free energy reaches its mini- 
mum (fic > M/j by increasing the current in the loop. 
Another important aspect to consider is that the three 
special magnetic moments, defined by Eqs. d^ . ([50]) and 
pip, are temperature independent. This means that the 
distance between two consecutive L and L + 1 magnetic 
moments, either free energy minima or zeros of the cur- 
rent, are temperature independent and equal to Ii /I2 and 
1 / Ii , respectively. For a very short ring, defined to have 
height comparable to the coherence length, Zq < ^0, the 
integrals Ii and I2 are easily computed, since it is enough 
to take the integrands at the center of the ring, z « 0, 
to obtain that /i ~ 77(0), and I2 ~ ?7(0)^. The three 
conditions, namely, the minimum of the free energy, the 
vanishing of the current and the highest LP tempera- 
ture take place at the same point /i//io = L/Ii because 



Ii = l2- Inserting the integrals and using the magnetic 
flux definition of Eq.(I71) in the center plane of the ring 
we obtain for the free energy density, critical temperature 
and current density the following expressions: 



\R 



I{l^,T) 



\T - 



Tc 
L - 

Tc{l^) 



<f(A.,0) 



$0 



Tc 



(f>(M,0) 



$0 



(33) 
(34) 

(35) 
(36) 



where we have defined the parameters, associated to the 
free energy, Tq = oi^/P, and to the current, /q, previ- 
ously defined. Notice that the above LP temperature 
was previously obtained using the linear approach. The 
above case is the standard LP effect, valid for a constant 
applied external field and arbitrary height. 

To show how the height of the cylinder affects the LP 
oscillations we discuss below the approximation that / 
is still a constant, but the integrals g and g' must be 
integrated along the z axis. We get analytical insight into 
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this problem by taking that the current loop is very small 
as compared to the thin-walled cylinder, a/R « 1, such 
that just the first term of the series expansion given in 
Eq.Q is enough to describe the external magnetic field: 



r]{z) = 



1 



El") 13/2 



(37) 



This vector potential is introduced into Eas. p3| and ([T 
to obtain that, 



(■ I 20-21 \ / ZO+21 \ 



} (38) 



2q — 21 



^ 20 + 21 - 



' 20-21 



20-21 ' 

R , 



, 2x 2 



' 20 + 21 ' 

V R . 



•J r 

+ — [ arctan 



arctan ( ^°^^' )]} 



(39) 



These integrals are expanded in a Taylor series in powers 
of [zqI , which is a small parameter, since the approx- 
imation of / constant is a good one only for a short cylin- 
der. We obtain for the integrals Ji and I2 the following 
values by keeping contributions of lowest order (zq/R)'^ 
and neglecting higher ones. 



Co 



1 



{1 



1 l-4(t)' 



2[l+(t)] ^ 



(40) 



We notice that the above equations do satisfy the in- 
equality of Eq.dUl), 



I2 -It ^3 



^o^2('2a^2/2l^2 
2 _ o\rJ \ rJ \ rJ 



2n5 



(41) 



As an example let us consider the short cylinder with 
radius R — lO.O^o, half-height zq — l.Ofo and the cur- 
rent loop just above its top, zi = l.S^o- Then h — Ii ^ 
( 10"^- The distance between two con- 

secutive free energy minima is 10.2456/io and between 
two consecutive current zeros is 10. 3507^0- 



IV. NUMERICAL ANALYSIS 

We obtain numerical solutions of the non-linear GL 
equation below, 

[ - + V{z)] f{z) = -a{T)f{z) -f Pf{zfm 



by solving it through the finite element method (COM- 
SOL and MATLAB softwares). Then we obtain the free 
energy, the current density and the current from Eqs.(IT|), 
([TT|) . and p^ . respectively, by input of f{z), obtained 
from Ea.(|42l). 

To obtain the numerical solutions of Eq. p2)) we cast 
it in dimensionless quantities, represented by a bar on 
top of the corresponding variable: / = •y//3/a(T)/, z = 



z/^T), where ^(r) = Co/Vl - T/T,, and V = V/aiT). 
The explicit temperature dependence disappears from 
Eq. (|42|) which becomes. 



(43) 



In practical grounds this means that all the lengths are 
measured in units of C(T), the temperature dependent 
coherence length, and so, ratios between two lengths be- 
come temperature independent. For instance, the func- 
tion ry, defined by Eq.(l5]), is found to be temperature 
independent, 77(2) = ?7(z), since it only depends on ra- 
tios. The dimensionless version of the potential, defined 
in Eq. (|44p . becomes temperature dependent and given 
by, 



V{-z) 



(44) 



where R{T) = R/^{T), fl{T) = m/Mo(T) and /io(r) = 
'i>o^{T)/2TT. Notice that the product no{T)R{T) = 
^qR/2tt is temperature independent. Like the radius, 
the height of the cylinder also shrinks to zero as T ap- 
proaches Tc'. zq = zo/£,{T). 

The free energy, the current density and the current, 
obtained from Eqs.dl]), ((TT|) . and ([12]), respectively, by 
input 01/(2:), determined from Eq. (j42|) . are also cast into 
dimensionless units. The current is given by. 



dz 



f{-zf[L--^{T) 



viz) . 
i?(T)J 



(45) 



where / = ///o(r), /o(T) = {qhZowa{T)) / {m(3R) . The 
free energy in reduced units becomes. 



F = 



dz rrdf(z) 

-20 ^0 



dz 



+ [v{z)-i]f{zr + -H-zr} 

(46) 



where F = F/Fo{T), Fo{T) = a(r)V2/3. 

We concentrate our numerical study in two thin walled 
cylinders, both with the same radius R — lO.O^o in pres- 
ence of the same current loop ring with radius a = 5.0^0, 
set at 0.5^0 above the top edge. The two cylinders 
only differ by the height, taken to be Zq = 2.0^0 and 
Zq — 10.0^0, such that zi = l.S^o and 5.5^0, respec- 
tively. We apply the same numerical method to both, 
but the short cylinder is expected to have the amplitude 
/ nearly z independent and be approximately described 
by the analytical ring limit previously considered. From 
the other side the tall cylinder, Zq = lO.O^Oj can only be 
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described by the full numerical treatment since the am- 
plitude varies along the cylinder's height, f{z). All the 
free energy and the current plots are in terms of the T = 
magnetic moment /i//io(0), in order to have a tempera- 
ture independent scale. Notice that fio{T) — ^q^{T) /2tt 
and we are taking its zero temperature value as our mag- 
netic moment scale. 

Figs. ^ and @ show the LP temperature versus mag- 
netic moment temperature diagram derived from the lin- 
ear theory. The ratio Tc(/i) /Tc is calculated from the low- 
est eigenvalue of Ea. (|i5| and for given L and fi/ /io(0) 
values. These L lines set the border line that separates 
the superconducting to the normal state in the diagram. 
The two cylinders display very distinct Tc(/i)/Tc versus 
^/^o(0) curves, as shown in Figs. ^ and 

Fig. ([2]) covers the to 200 /i//io(0) range, which con- 
tains the first 20 L lines of the Zq = 2^o cylinder. Notice 
the parabolic shape of the L lines and their nearly ho- 
mogeneous spatial distribution. The distance between 
two consecutive maxima, of L and L + 1 lines, is de- 
scribed by the formula A/i//io(0) ~ 9.5. Notice that for 
increasing magnetic moment the ratio Tc{p)/Tc becomes 
lower than one, indicating that the height of the cylinder 
matters in this case, and our theoretical considerations 
developed in the ring limit applies to this situation and 
is determined by Eg.ipUl and not by Eq. ([Ml) . The ratio 
Tc{ii)/Tc reaches one only in case there is only a single 
magnetic flux piercing the ring, that is, in case of the 
nearly zero height cylinder or of the constant magnetic 
field applied to cylinder of a arbitrary height. 

Fig. [3] shows a much richer structure for the L lines 
of the R — 10^0 cylinder, due to the inhomogeneity of 
the magnetic field that results in distinct magnetic fiuxes 
piercing this cylinder, each for a different z plane inter- 
sect. Notice the presence of an inner and an external 
(true) border lines. To understand the R = lO^o dia- 
gram, we selected a few L lines, plotted as thick lines 
for best visualization purposes: L = 5, 12, 20 and 41. 
The L — 12 line is a key case because of its double peak 
structure with a local minimum between them, a com- 
mon feature to many lines. The first peak is narrow and 
tall, while the second peak is short and broad and takes 
place at a larger magnetic moment than the first peak. 
The tips of the narrow and broad peaks of the L = 20 
line still fall inside the plotted temperature range, 0.8 to 
1.0 Tc{iJ.)/Tc. The L = 5 thick line shows the onset of 
the second peak, not present in the lower L lines, which 
only have a single peak. The L — 5 and 6 can be re- 
garded as transition lines to the onset of the second peak 
and L = 7 is the critical line because the second peak is 
clearly established for L = 8 and beyond. L = 41 is the 
highest line in this diagram, seen around /i//j,o(0) = 400 
and setting the end of an inner transition line inside the 
superconducting region. The first (narrow) peak con- 
tributes to the true border line only up to L = 4 then 
to submerge and form an inner border line for high L 
values. Thus this inner border line is the curve tangent 
to the maximum of the first peaks. From its side, the 



second (broad) peak only contributes to the true super- 
conducting normal state border line, from the L — 5 line 
up to i = 20 line. For the sake of the argument we have 
included the L = 5, 6 and 7 lines in this set, although 
the second peak is only in an embryonic level there. In 
conclusion the superconducting normal state border line 
has less angular momentum states for the tall cylinder, 
since L = 20 is the maximum possible value reached at 
/i//io = 510, whereas for the short cylinder, L = 20 is 
reached for fi/fio — 200. 

The full non-linear theory is used to obtain the next 
plots, which means that Eq. (j43p is numerically solved by 
means of the finite element method, and next the free 
energy and current calculated from Eqs. (H51) and (^5]) . 
respectively. In the dimcnsionless treatment the temper- 
ature T is introduced by adjustment of the radius and of 
the height, given by R{T) and Zo{T), respectively. 

Fig. m shows several L free energy lines of the Zq ^ 
2-0^0 cylinder for the two temperatures of T/Tc = 0.8 
and 0.9. The minimum of each L line is also a free en- 
ergy minimum. These minima are equally spaced and 
two consecutive ones, associated to L and L-\-l lines, are 
separated by A/i//io(0) ~ 9.5 in agreement with our the- 
oretical analysis that predicts A^/^o(0) = Ii/h, accord- 
ing to Eq. ip^ . This position is temperature independent 
and Fig. |4] has a dashed straight line to indicate that for 
the two temperatures the minimum is in the same mag- 
netic moment value. 

Fig. [5] shows the vanishing of the current of each L 
line for the two temperatures of T/Tc = 0.8 and 0.9. 
The zeros are found to be equally spaced, A/i//io(0) ~ 
9.5, also in agreement with our theoretical analysis that 
predicts A/x//xo(0) = l//i, according to Ea. (pO)) . The 
dashed vertical straight line is a guide to the eye to show 
that the current vanishes at the same magnetic moment 
independently of the temperature. 

Fig. ini shows the free energy of the tall cylinder {Zq = 
lO.O^o) for the temperature of T/Tc = 0.8. Distinctly 
from the short cylinder case, here the minimum of a i 
line is not a free energy minimum, unless approximately 
for the first few lines, L = to 4. The L lines undergo a 
transition in shape from the small to the large L limits. 
To understand this transition we selected the L — 3 and 
12 lines, also shown in the inset, as they are represen- 
tative of each side of the transition. While the L ~ 3 
line has a single minimum, L = 12 has double well shape 
with two minima at fj,/fj,o{0) ~ 130 and 270, the first 
one being the lowest one in free energy. Nevertheless the 
second well is more important than the first one because 
it sets the true free energy minimum at /i//io(0) « 300, 
thus away from the L line secondary local minimum. The 
first one has other competing states with lower energy, 
the lowest one being the L = 7 state, which is drawn as 
a thick line to help the visualization. The L = 7 state is 
a critical one as it sets the onset of the double well struc- 
ture. This double well structure is not a feature present 
in all lines, the first ones simply do not have it. In conclu- 
sion the free energy clearly shows two types, single and 
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double well shaped, that dominate the small and large 
magnetic moment regimes, respectively. This behavior 
is also seen in the temperature T/Tc = 0.9, as shown in 
Fig 13 There L = 5 is the critical line, and we have also 
plotted as thick lines, the L — 3 and L = 8, to show the 
single to double well transition. 

Fig. [8] shows that the Zq = lO.O^o cylinder displays 
novel features as compared to the standard LP problem. 
The current of the transition L line that separates the 
single (low L) to the double well (high L) behavior dis- 
plays unusual behavior. Both limits are characterized by 
distinct slopes but apart from this both single and dou- 
ble well regions are similar, as shown for the , but not in 
the transition (L = 7) line. This is seen for both temper- 
atures, T/Tc — 0.8 and 0.9. The transition lines occur 
for L = 7 and 5, respectively and displays a very unusual 
current pattern near its zero. 

The zeros of the current and the free energy minima 
coincide for the standard LP, but the inhomogeneous ex- 
ternal field sets them apart, as shown here. For this 
reason we plot the free energy for a given temperature 
(T = 0.8Tc) together with straight vertical lines associ- 
ated to the zeros of the current. Fig. [S] shows for the 
Zq = 2.0^0 cylinder that these straight vertical lines fall 
very close to the free energy minima in agreement with 
our theoretical predictions derived in the ring limit. We 
have selected a few magnetic moment values (red dashed 
vertical lines) in this diagram for further analysis, given 
by fi/fio = 9.62, 66.64, 114.04 and 190.29, whose free en- 
ergy minima correspond to the L = 1, 7 12 and 20 lines. 
Interestingly we find that these vertical lines intercept the 
matching points between the L — 1 and L + 1 free energy 
lines up to numerical precision. The distance between 
consecutive vertical lines is in agreement with the theo- 
retical prediction for A^c/mo(0) « 9.5. We observe that 
this property seems to be general, namely, also valid for 
the Zq = 10.0^0 cylinder, as shown in Fig. (TU] However 
for the tall cylinder the situation is quite different as the 
vertical lines clearly form two sets, a low and a high mag- 
netic moment set with a gap between them, located at the 
L = 7 line. The two sets feature distinct separation be- 
tween consecutive vertical lines with very distinct behav- 
ior regarding the free energy minima. The vertical lines 
of Fig.fTUlsplit into two sets according to their separation. 
For low values of magnetic momcnt(until L = 5) the dis- 
tance between two consecutive lines is A/i//j,o(0) « 14.7. 
Then, there is a transition region(L — 6,7,8), and for 
L > 9 the curves of the free energy are dominated by 
the second well. In this region, the distance between the 
lines is given by A/Lt//Lto(0) ~ 30.1. Comparison betwee 
tall and short cylinders show that their low magnetic mo- 
ment has similar behavior, but not for the high magnetic 
moment regime because of the double well shape of the L 
lines for the tall cylinder. The zeros of the current are far 
from the free energy minima and yet they coincide with 
the crossing of the L—1 and L+1 free energy lines. In this 
case we have also selected a few magnetic moment val- 
ues (red dashed lines) for further analysis: fJ,/fJ,o = 12.75, 



85.64, 129.02 and 264.97, They are labeled 1, 6, 7 and 11 
because they fall within the regimes where these L lines 
provide the free energy minima. 

In the ring limit (very short cylinder) it is easy to verify 
that the crossing point between any L ~ I and L + 1 free 
energy curves happens at the magnetic moment value of 
the vanishing of the L state current. In the ring limit the 
amplitude / is constant and the free energies Tl-i and 
J-L+i for states L—1 and L+1, respectively, are given 
by. 



-^[«o(|)V-i+«(r)]^ and (47) 
-^[«o(|)^.9L+l+«(^)]^ 



(48) 



according to Eq. p^ . Next take that J-l-i — to 
obtain that gL~i{p*) = gL+i{^*), which gives that. 



zO 



^0 



L + 1 



L - 1 



$0 



$(^*,z) 



that once simplified becomes. 
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zO 



dz 
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$(/x*,z) 
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(49) 



(50) 



Inserting the definition of the magnetic flux $(/^*,z) 
(Eq.Q), and expressing the result as a function of in- 
tegral /i, previously defined, we obtain that. 



(51) 



This shows that for constant amplitude / it holds that 
^* = A^c, according to Eq. (pO)) . We also stress that our 
present numerical evidence shows that this property also 
holds for the tall cylinder where the amplitude / is no 
longer constant. 

To gain further insight into the structure of the su- 
perconducting state along the cylinder height, we plot in 
reduced units the current density J{z), the superconduct- 
ing density |-0(z)p, and the potential V{z) as a function 
of z. Notice that the current loop, which is the source 
of the inhomogeneous magnetic field, is located in the 
positive z axis 0.5^o above the maximum height. Fig. 1111 
shows these plots for the Zq = 2.0^o cylinder at the se- 
lected magnetic moment values called (1), (7), (12) and 
(20). All the possible crossing L states with the verti- 
cal lines, as shown in Fig.|9l are in the J(z), |V'(^)P and 
V{z) plots. To understand the information given by these 
plots we start the discussion of the J{z) plots. Because 
these magnetic moment values are very near to the zeros 
of the current we take this property to claim that the 
sought free energy minimum must be such that its J{z) 
integrated along the z-axis will add to a zero current. 
Fig. ITlT la') shows that only the L = 1 curve qualifies 
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for a free energy minimum since it is the only one to 
cross the z-axis thus yielding positive current for z < 
and negative for z > 0. According to this criterion all 
others L lines can be dismissed as they can not be zero 
current states since they do not cross the z-axis. This 
reasoning also provides the key to select the free energy 
minimum in Figs. Illf 2a). (3a) and (4a). The present 
criterion uniquely selects the L = 7, 12 and 20 states 
among the other ones drawn in these figures, respectively, 
as the only possible free energy minimum because only 
them fit as zero current states. For the short cylinder the 
superconducting density |'0(z)p, and the potential V{z) 
also provide features able to uniquely define the sought L 
lines. For the free energy minimum states the density is 
nearly constant according to Figs. ITlT lb). (2b), (3b) and 
(4b). Similarly the potential approaches zero for these 
states as seen in Figs. ITlT lc'). (2c), (3c) and (4c). 

FigdU shows these plots for the = lO.O^o cylinder 
at the selected magnetic moment values called (1), (6), 
(7) and (11). Fig. [T^ la) shows that in the low mag- 
netic moment regime the L state associated to the free 
energy minimum can be judiciously selected by the zero 
current condition and this is the L = 1 state. The same 
does not hold for the intermediate magnetic moment val- 
ues (6) and (7) where many lines qualify at least at the 
naked eye level precision, which means that many lines 
cross the z-axis near the center of the cylinder having 
positive current for z < and negative for z > 0. For 
the high magnetic moment limit the situation is simi- 
lar although the current vanishes near the upper edge 
of the cylinder, z ~ 5, because of the proximity to the 
strong external current loop. The dramatic changes that 
the superconducting state undergoes by increasing the 
magnetic moment are witnessed in the superconducting 
density |?/'(z)P, and the potential V{z). Figs. HUlb), 
(2b), (3b) and (4b) show that density is mostly concen- 
trated at the top edge of the cylinder for fx/ = 12.75 
which falls in the low magnetic moment regime. For 
the intermediate values of /i//io = 85.64 and 129.02 we 
find L states concentrated at both edges to finally find 
that in the high magnetic moment regime, represented 
by /i//io — 264.97, the density is mostly concentrated at 
the bottom. Indeed Figs.lT^lc). (2c), (3c) and (4c) show 



that the potential undergoes dramatic changes from the 
low to the high magnetic moment regime. In the former 
case (/i//io — 85.64) it undergoes mild changes from one 
edge to the other as it varies from to 1 in reduced units, 
while for the latter case(/i//io — 264.97) it varies from 
to 23, being very intense in the top edge. These plots ex- 
plain the reasons for the twofold regimes present in the 
LP oscillation near the persistent current loop. 

V. CONCLUSION 

We have considered the Little-Parks oscillations on a 
thin-walled superconducting cylinder with a persistent 
(magnetic moment) current loop set on its top. This 
mesoscopic system is shown here to have two regimes of 
temperature oscillations not found in the original Little- 
Parks system because instead of a constant applied mag- 
netic field there is an inhomogeneous magnetic field cre- 
ated in space by the persistent current loop. We ob- 
serve that this inhomogeneity grows in relevance accord- 
ing to the cylinder height, as expected. In order to un- 
veil the new features we have considered a short and a 
tall cylinder, and report new results in both cases. For 
the short cylinder we report new analytical expressions 
for the magnetic moments associated to the free energy 
minima, to the zeros of the current, and to the maximum 
attainable temperature (Eqs.d^Hl), (EOl) and ^Q). For in- 
stance, we find that the current in the thin-walled cylin- 
der vanishes at a magnetic moment greater than that of 
the free energy minimum (/ic > M/)- For the tall cylinder 
we find novel features, perhaps best captured by Fig. [S] 
This figure shows two distinct current regimes for the tall 
cylinder, set at low and high L, respectively, associated 
to distinct slopes with a clear transition L line that sep- 
arates them. The existence of these two regimes, and the 
transition that separates them, are within the realm of 
observable experimental measurements. 
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